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The space radiation problem 



□ 


□ 


□ 


□ 


Space radiation is comprised 
of high-energy protons and 
heavy ions (HZE’s) and 
secondary protons, neutrons, 
and heavy ions produced in 
shielding 

Unique damage to bio- 
molecules, cells, and tissues 
occurs from HZE ions that is 
qualitatively distinct from X- 
rays and gamma-rays on Earth 

No human data to estimate risk 
from heavy ions, thus requiring 
use of biological models and 
theoretical understanding to 
assess and mitigate risks 

Shielding has excessive costs 
and will not eliminate galactic 
cosmic rays (GCR) 



iron 
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Cucinotta and Durante, Lancet Oncology (2006) 




Radiation tracks and energy deposition 

□ The energy deposition by heavy ions is highly heterogeneous and 
dependent on the type and energy of the ion 

□ The interactions of radiation with matter are stochastic in nature and 
therefore often studied by Monte-Carlo simulations 

Primary energy loss events in Primary energy loss events in 
low-LET tracks high-LET tracks 




~ioo A 

ClOOeV 


-100 to 500 eV 
blobs 


Short tracks 

-500 eV to 5 
keV 


Primary 

> 


Branch tracks 

>5 keV 


Delta rays 


A. Mozumder and J.L. Magee (1966) Radiat. Res. 28, 203 


Penumbra 


Track core 


Delta rays 

C. Ferradini (1979)/. Chim. Phys. 76, 636 
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Radiation effects: time sequence of events 



Time (s) 


Stage 


Events 


Modeling 


io - 15 

10-12 


Physical 

Physico- 

chemical 


Energy absorption 

Reorganization 
Electron thermalization 


Particle transport 
Cross sections 


10-9 

10-6 

io - 3 


Chemical Radical diffusion Green’s functions 

Chemical reactions 

Biological DNA repair Kinetics models 

Molecular dynamics 


V 
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Radiation effects: time sequence of events 
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Particle transport basics 


□ The trajectory of a particle and all its interactions is 
followed in the medium 

□ A particle is followed until 

• It leaves the volume of interest 

• Its energy decrease below a threshold 

• It disappears by a physical process (e.g. absorption of a photon 
during a photo-electric effect 

□ Many other particles are generated by the interactions 
of the “primary” particle. The trajectories of these 
secondary particles should also be followed. 


Page No. 6 Plante, I. and Cucinotta, F.A. (2011) Monte-Carlo Simulation of Ionizing Radiation Tracks. In: Mode, C.J. (Ed.) Applications of 
Monte Carlo Methods in Biology, Medicine and Other Fields of Science. InTech, Rijeka, Croatia, www.intechopen.com 


Particle transport basics 


v = 


sin(O) cos(cp) 
sin(0) sin(cp) 
cos(O) 


□ Particles 

• Position (x,y,z) 

• Energy (E) 

• Direction (0,cp) 

□ Cross sections 

• Probability of interaction between radiation and matter 



dl = -Incdx 


I: Incident fluence 
n: Density of targets 
dx: Width 
a: Cross section 


Incident 

fluence 


dx 



rv . 


0 \ 


© ® T 



0 

© 



l® © 


a ^ 

r 


Transmitted 

fluence 


Cross sections (units: cm 2 ) 
Mean free path 1 (units: cm) 


I(x) = 1(0) exp(-x / 1(E)) 
1(E) = 1 /(Ng(E)) 
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Particle transport basics 


□ Cross sections (total and differential in energy, angle,... i.e. do/dW, 
dcj/d0, d 2 a/dWda,...) are needed for particle transport 

□ The total cross section is the sum of the cross sections for each 
interactions 

□ A distance to the next interaction (s) is sampled from the exponential 
distribution: 

s = —X log(U) 

ME) = 1/(Ng(E)) 

□ The particle is moved to the next interaction site: 

x'= x + ssin(0)cos((p) 

y’ = y + s sin(0) sin(cp) 
z'= z + scos(0) 

□ The phenomenon (e.g. ionization) is determined by the ratio of the 
cross section a phenomenon over the total cross section 

□ The energy loss and change of direction are determined by sampling 
the differential cross sections 
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Cross sections 


□ RITRACKS includes accurate cross section models for 
all ions and secondary electrons or photons 


□ For electrons: 

• Ionization 

• Excitation 

• Elastic collisions 

• Dissociative electron 
attachment 

• Bremsstrahlung 


□ For ions: 

• Ionization 

• Excitation 

□ For photons: 

• Compton effect 

• Coherent diffusion 

• Photoelectric effect 

• Pair production 
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Cross sections 


Cross sections 


□ Cross sections used in RITRACKS 

Ions 



The cross sections for ions are scaled with Z 


eff- 


do io„ ( V ) _ <72 da prolo„( V ) 

” ^eff 


dW 


dW 


Z eff /Z = l-exp(-125p 2 /Z 2 ' 3 ) 


Electrons 



v: velocity of the ion 
(3: relativistic v/c 
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Monte-Carlo simulations 


□ Used to simulate stochastic systems 


N=10000; N = 7863; ti=4*N /N=3.1452 




=>By increasing the number 
of random points, 4N C /N 
converges to the value of n 


Example: Calculation of n 
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Monte-Carlo simulations 


• Random numbers 

• Probability distributions 

• Normalization 


N 

I>i =1 

i=l 

00 

Jf(x')dx'= 1 

—oo 


Discrete 


Continuous 


Uniform distribution 



Poisson distribution 


P(x=k)=A, x e'Vx! 



Exponential distribution 
p(x)=Ie' Xx 
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From http://en.wikipedia.org/wiki/List_of_probability_distributions 


Monte-Carlo simulations 


□ Generating probability distributions 




Example: Exponential distribution 

f(x)=^e'A A>0, x>0 
F(x)=l-e^ x =U 

X=-Log(l-U)A or X=-Log(U)A 

U is a random number between 0 and 1 
F(x) is the cumulative probability distribution 

X 

F(x) = |f(x’)dx' 

k 

P k =IPi 

i=l 
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Monte-Carlo simulations 


U 


□ Generating the exponential distribution 


Sampling of the exponential distribution 
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The software RITRACKS 



□ The software RITRACKS comprises several parts 

• The calculation part includes: 

S The cross sections, which are necessary for particle transport 
S The particle transport routines 
S Post-simulation data management 

• The Graphic User Interface (GUI), comprises several windows: 

S The main window 

S Incident radiation window 

S Cross sections windows (electrons and ions) 

S Results (events) details 

• The 3D visualization window 

• The help file 

□ All necessary files are included in an installer for 
Windows 
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The software RITRACKS 



□ Multiple CPU computing 
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The software RITRACKS 


□ Multiple CPU computing (Windows) 


RITRACKS VI 


RITRACKS V2 
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Heavy ion track structure simulation 


Simulation fo 1 H + , 12 C 6+ , 28 Si 14+ and 56 Fe 26+ tracks, 100 MeV/amu 
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Dose (Gy) Dose (Gy) 


Radial dosimetry 




(a) 'H + 1 MeV/amu, LET-33 
keV/ F m 

(b) 20 Ne 10+ 377 MeV/amu, 
LET-31 keV/nm 

(c) 4 He 2+ 1 MeV/amu, 56 Fe 26+ 1 
GeV/amu, LET-150 keV/ F m 
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Voxel dosimetry 



1 GeV/amu 56 Fe 26+ ion 

LET-150 keV/nm 

Voxels: 40 nm x 40 nm x 40 nm 
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Target dosimetry 




X (nm) 


Page No. 21 
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Target dosimetry 
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DNA damage / yH2AX foci studies 


• Irradiation by 1 
GeV/amu Fe ions 

• 1 00 cGy 

• LET - 149keV/|jm 


Tracks only 


Nuclei only 


Tracks and nuclei 


'! A- 
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Experiments performed at the NASA Space Radiation Laboratory (2007) 







DNA damage 


Calculated DSB 


/ 




2700 x 1 H + , 300 MeV (1 Gy) 
Dose in voxels (20 nm) 
Chromosomes (RW model) 
Intersection voxels 
H2AX foci experiments 
Application of DSB probability 
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DNA damage 


Calculated DSB 



6 x 56 Fe 26+ , 1 GeV/u (1 Gy) 
Dose in voxels (20 nm) 
Chromosomes (RW model) 
Intersection voxels 
H2AX foci experiments 
Application of DSB probability 
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DNA damage / yH2AX foci studies 


• Calculation of DSBs by low- and high- 
LET radiation 
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DSB/bp 


DNA damage / yH2AX foci studies 


Calculation of DSBs by 1 H + , 12 C 6+ and 56 Fe 26+ 
ions 


Protons 1 GeV 


Carbon 293 MeV/u 


Iron 1 GeV/amu 
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DSB/Gy/bp 


DNA damage / yH2AX foci studies 


• Calculation of DSBs vs LET by 1 H + , 12 C 6+ and 
56 Fe 26+ ions 



3 


LET (keV/|im) 


LET (keV/nm) 


LET (keV/nm) 
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ti DSEs per II 


DNA damage / yH2AX foci studies 



• DSBs per cell vs dose is linear 

• Exchanges per cell vs dose is linear-quadratic 
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From Dr. Artem Ponomarev 


DNA damage simulations 



□ To better understand the formation of DSBs, a chromatin 
fiber is build from nucleosome units and linker DNA 

Nucleosome (DNA fragments) Chromatin fiber Nuclei 





Nuclei simulations courtesy of Dr. Artem Ponomarev 
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DNA damage simulations 


□ The Binary-Encounter-Bethe (BEB) model of ionization cross section 


□ The energies are expressed in units of ionization potential of the 
orbital (B): 

• t=T/B is the kinetic energy of the incident electron 

• w=W/B is the kinetic energy of the ejected electron 

• u=U/B is the kinetic energy of the electron in the orbital 

□ The total cross section is obtained by integration 


dw t + u + l ( t — w ) 2 (1 + w) 


da s 1 1 


l 


1 

l + t t — w 1 + w 



a = 
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DNA damage simulations 



□ In the DNA bases, there are many internal and valence electrons. The BEB 
model allows to model the ionization for each electron of the molecule. 


Thymine 



No 

U (eV) 

B (eV) 

N 


1 

794.1 

559.41 

2 

00 

2 

794.1 

559.18 

2 


3 

601.5 

425.6 

2 

c 

o 

4 

601.5 

425.48 

2 

u 

_QJ 

5 

435.9 

311.36 

2 


6 

435.9 

310.43 

2 

c 

cu 

7 

435.8 

308 

2 

c 

8 

435.7 

306.41 

2 


9 

435.8 

305.8 

2 


10 

66.26 

40.06 

2 


11 

71.74 

39.14 

2.04 


12 

63.05 

36.16 

1.86 


13 

57.89 

34.56 

1.85 


14 

44.95 

30 

2.3 


15 

43.96 

26.22 

2.35 


16 

47.64 

25.14 

2.11 


17 

40.64 

24.85 

2.15 

'T 

18 

41.92 

21.32 

2.18 

o 

19 

39.28 

20.94 

2.02 

U 

20 

36.32 

19.4 

1.86 

QJ 
C V 

21 

44.53 

18.7 

2.05 

<1> 

U 

C 

22 

55.89 

18.56 

1.92 

75 

> 

23 

54.72 

17.59 

1.99 


24 

39.88 

16.62 

2.01 


25 

46.93 

16.15 

2.03 


26 

39.89 

15.44 

1.78 


27 

47.3 

13.96 

1.83 


28 

59.96 

13.15 

1.86 


29 

54.12 

12.31 

2.07 


30 

60.23 

12.10 

1.78 


31 

40.38 

9.27 

1.97 



thymine. log 
Model 1.2 _ MO 1/147 
Energy = -19.L29S2 




thymine. log 

Model 1.2 _ MO 10/ 14 7 

Energy = -1.07502 
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DNA damage simulations 


□ The BEB cross section can be written in a form suitable for sampling 
by a composition method: 


da 

dw 

da 

dw 


+ 


t + u + 1 |^(Y — w)^ (1 + 

Z k i = Z A tSi ( w ) h i o) 


1 + t 


1 1 
+ 


t — w 1 + w 


+ 


/=! 


/=! 


1 


+ 


1 


(t - w) (1 + w) 


log(0 


where 



4 s 

' 1 


( 7 - 1 ) / 2 

^1 _ 

£ + W + 1 

+ lj 

+iy 


jg,.(w)Jw = l 
0 

0 < h i (w) < 1 
A : >0 


i 

r ^ ^ 

/-I 

2 

^ t + + 1 y 

t(t + 1) 


K (w) = 


t + 1 


w + 1 
t + 1 


h 2 (w) = 2 


1 - 


t-w 
( + 1 


g i(w) = 


g 2 ( w ) 


t + 1 1 

r-iO+i) 2 

((1 + 1 ) 1 
t — l ( t — wy 


A 3 =2 log (0 


f ^ 1 

"0 + l) 2 -4" 

A (w) = — 

1 + 

f w + 1^ 

3 

g 3 ( w ) = 

<N 

? ; 
cT 

1 

^ + W + 1 y 

1 

(N 

+ 

cnT 
1 

3V 7 2 


[t-wj 



L^+ir-4j 


1 


(w + 1) 3 
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DNA damage simulations 


□ Sampling algorithm and results 

Algorithm to sample the energy of the ionized electron w* from the BEB model of ionization 
cross sections 


CALCULATE A^, A 21 A; and A^25 — A1+A2+ A3 
REPEAT-; 

GENERATE random numbers R lf P 2 and R = ; 

IF (R 1 <A 1 /A 1Z 3j { 

{ 

. w-* 

' 9 + i)-«,(t-i) 

t — -ft! 

IF ( R-i < ) EXIT loop 


} 

IF A 1/ ^ 1 2 i, Ai+A;)/ Aj_ 23) { 


w* = 


' (f+l) + A(f-l) 


IF [R 3 <2 


' l -till' 

t + 1 


1 EXIT loop 


} 

IF (Ri>(Ai+A 2 )/A 1h ] 


t + 1 


W* = 


IF i A < 


> + CA-i)Ca + i) 2 -4) 

1 


-1 


” 

f ! \ 

r~ 

TO 

1 + 

w, + l 



~ J 

1 


EXIT loop 


} 

RETURN 


Sampling of the BEB differential cross sections of Thymine 
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DNA damage studies 



□ Cross sections can be calculated for the bases, sugars and phosphates. 

□ In this case, the medium is considered a succession of homogeneous 


media. 

dl = -INadx 


X 

ln(7 / 7 0 ) = - f <s{u)Ndu 


I - I 0 exp 


N 


o 

/— i 


ZK -OiWj-KSiX 

7=1 


Relative weight 


z-1 


exp(-7V Z WjGj )(1 - Qxp(-NW i o i )) 

D = I 1 — 

Pi 1 0 


Mr. 


Sampling of W s 


w s = 2 Wj - -±- log[l - V(1 - e °‘ NW ‘ )] 


7=1 


No, 


Oi 


-> 

-> 


> 


-> 


G c 


Gc 


g 4 


< — W 1 — > 


W 2 -H 


<- W 3 -> 


«- 


w. 


Relative intensity 



■> 
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Radiation chemistry 


□ ~10- 12 - 10- 6 s 

• Particles diffusion 

• Chemical reactions 


Number of chemical species created 

G(X) = 

100 eV deposited energy 


□ The radiolytic species are not uniformly distributed. 
Therefore, an approach based on Green’s functions of 
the diffusion equation (DE) is used. 


Examples of chemical reactions: 


+ e- 2H 2° 

c aq c aq ► 

H 2 + 2 OH 

•OH + e-„ — * 

OH' 

’OH + ’OH * 

H 2 0 2 


’oh + h 2 o 2 — ► ho 2 * + h 2 o 

H + + 0 2 *' H0 2 * 

H’ + H 2 0 2 * ’OH + H 2 0 
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(More than 60 reactions...) 


Bimolecular reactions 


□ DE for the propagation of particles A and B 


5p(r A ,r B ,t|r A0 ,r B0 ,t 0 )/a 


d a v a 2 + d b v b 2 


]p( r A» r B> t I r AO> r BO’ t o) 


p(r A ,r B ,t| r Ao r Bo ,1^,): probability of 
particles A and B to be at positions r A 
and r B at time t, given that they were at 
r A0 and r B0 at ^ 


D a , D b : Diffusion coefficients 

□ Transformation of variables 


R - (D B r A + D A r B )/(D A + D b ) 

r = r B“ r A 


V R =d/dR 
V r = d/dr 


<3p(R>r,t | R 0 ,r 0 ,t 0 )/5t = (D a +D b )[v r 2 + V 2 |}(R,r,t | R 0 ,r 0 ,t 0 ) 

p(R,r,t | R 0 ,r 0 ,t 0 ) = p R (R,t | R 0 ,t 0 )p r (r,t |,r„,t 0 ) 


3p R (R, 1 1 R„t 0 )/5t = (D a +D B )V„ 2 p R (R,t | R„ t 0 ) 


Uncoupled equations 


0p r (r,t | r 0 ,t 0 )/a = (D A +D B )v 2 p r (r,t | r 0 ,t 0 ) 


in r and R 
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Bimolecular reactions 


□ Free diffusive motion of the coordinate R 


5p R (R,t | R o ,t o )/0t = DV R 2 p R (R,t | R 0 ,t 0 ) 


(DE) 


p R (R,t | R 0 ,t 0 ) -8(R-R 0 ) 


P R (I R |— > °°, 1 1 R 0 ,t 0 ) — 0 
p R (R,t|R 0 ,t 0 ) = 


1 


r A T'X / a , XT 3/2 


[4jrD(t-t 0 )] ; 


(R-Rg ) 2 

4D(t-t 0 ) 


(Initial condition) 
(Boundary condition) 


(Solution) 


p R (R,t|R 0 ,to): probability distribution of the vector R at time t, given that it 

was located at position R 0 at time t 0 

D=D a +D b : Sum of the diffusion coefficients 
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Bimolecular reactions 


□ Free diffusive motion* of the inter-particle separation vector r 


dp r (r,t|r„,t 0 ) 

St 


DV r 2 p r (r,t|r 0 ,t 0 ) 


(DE) 


P r (r,t|r # ,t 0 ) = 5(r-r # ) 


(Initial condition) 


□ For chemical reactions, we need the inter-particle distance r. 

• Therefore, the DE is written in spherical coordinates. 

• Only the radial component will be considered (angular dependency terms are neglected). 
This considerably simplifies the analytical solution. 


dp( r >t | r 0 ) 
St 


D_8_ 
r 2 Sr 


0 , , x 

— P(r, 1 1 r 0 ) 
or 


4 nr 0 2 p(r, 1 1 r 0 ) = 5(r - r 0 ), r > R 


(Radial part of the DE) 

(Initial condition; R=reaction radius) 


p(r,t| r 0 ,to): probability distribution of the separation distance r at time t, given that it was r 0 at time t^ 
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Bimolecular reactions 


□ Simple case: partially diffusion-controlled reaction 
with rate k a and reaction radius R 

k a 

A + B^C (Reaction) 


\nR D 


2r ,ap(r,t|r 0 ) 


dr 


= k a p(R,t |r 0 ) 


r=R 


(Boundary condition) 


47nr 0 p(r,t|r 0 ) = 


a / 47iDt 


i exp 


(r-r 0 ) : 

4Dt 


+ exp 


(r + r 0 -2R) 2 
4Dt 




+ aW 


V 


r + r 0 - 2R 
V4Dt 


,-aVDt 


(Green’s function) 


Q(t | r 0 ) = J47tr 2 p(r,t | r 0 )dr = 1 + 

R 


Rot + 1 
r 0 a 


W 


r o -R 

V4Dt 


5 a 


VDt 


Erfc 


t r 0 -R 

V4Dt. 


(Survival probability) 


k n + 47rRD 

a = — 

4ttR 2 D 

The probability of reaction P(t|r 0 ) = 1 - Q(t|r 0 ). At each time step, the probability of reaction is assessed. If 
the particles have not reacted, their relative distance is obtained by sampling the Green’s function using 
the algorithm described in Plante et al. (2013). 
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47ir 2 p(r,t|r 0 ) 


Green’s function for radiation chemistry 


□ Simple case: partially diffusion-controlled reaction with rate k a and 
reaction radius R 

k 

A + B -> C 


Green’s functions 



Survival probabilities 
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Green’s functions (3D) for radiation chemistry 


□ Green’s function with angular dependency (exact) 

p(r, 0, <M | r 0 , 0 O , <|> 0 ) = - — i= Yj ( 2n + (cos(y)) f e"" 2D 'uF n+1/2 (u, r)F n+1/2 (u, r 0 )du 

4n^rr 0 „ =0 o 

F ( . T (2RF +l)[l v (ur)Y v (uR)- Y u (ur)J u (uR)]-2uR[j v (ur)Y;(uR)- Y v (ur)JUuR)1 
V[(2Rk a + 1) J v (uR) - 2uRj;, (uR)] 2 + [(2Rk a + 1) Y v (uR) - 2uRY; (uR)f 


cos(x) = cos((9)cos((9 0 ) + sin((9) sin((9 0 ) cos(^ ~(/> 0 ) 


P n (x) n th Legendre polynomial 

J n (x) Bessel function of the first kind 
Y n (x) Bessel function of the second kind 

=> Even though this distribution is very complex, it is possible 
to generate the exact length and the approximate angles of the 
inter-particle vector from this Green’s function 
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Green’s functions (3D) for radiation chemistry 


• Limiting case: R— >0, k a =0 

• In principle, this should be equivalent to free diffusion 


p(r, e, <M I r 0 ,e 0 ,<|> 0 )= . \ — Z ( 2n + ! ) p n ( cos (t)) J e VDt uF n+ 1 n (», r) F n+ i , 2 (». r o )du 

A n-r / w* * 


47r A /rr n n =0 


F (1 , r) _ [Jy (ur)Y v (uR) - Y„ (ur) J„ (uR)]- 2uR[j v (ur)Y; (uR) - Y v (ur) J’ v (uR)] 

Vhv (uR) - 2uRj;, (uR)f + [y v (uR) - 2uRY; (uR)f 


lim F u (u,r) = J u (ur) 


0 


P (r, 0,<f>,t\v o ,0 o ,</> o ) = — — e - (r2 ^ ),4D ‘ 


4 TT-JrrQ 2Dt 


£(2w + l).P„(cos(;K))/, 


IT„ ' 


n+ 1/ 2 


n = 0 


\2Dt j 


Weber’s formula 


j^'Y( 2 k + 1 )I t+1/2 (z)P k (cos y) = e 


zcosy 


p(r,0,</>,t\r o ,0 o ,0 o ) = 


1 


-(r 2 +r 0 2 -2rr 0 cos /)/ 4 Dt 


(4 7rDt) 


3/2 


= Free diffusion! 
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Green’s functions (3D) for radiation chemistry 


Radial distribution of the exact Green’s function 

The radial distribution is not function of 0 O and (|) 0 , so we use 
0 o =O and 4* 0 = 0, for which cos(y)=cos(0) 

In 7i n 

p(r, 1 1 r 0 ,0,0) =|| r 2 p(r, 6, </>, 1 1 r 0 ,0,0)sin6ti6ti^ = J 2^r 2 p(r, r 0 ,0,0)sin6ti6 > 

0 0 0 


n 

J (2« + 1 )P n (oos6)s\n(0)dO = 2 S nC 


Identity 


2 00 


p(r, 1 1 r 0 ,0,0) = -= £ S n0 J e' u2Dt uF n+1/2 (u, r)F n+1/2 (u, r 0 )du = ~^= J e‘ u2Dt uF 1/2 (u, r)F 1/2 (u, r 0 )du 

V^o »=0 0 


0 0 


J„„ 2 (z) = V j 0 (z) = 


sin z 


Y n+ 1 / 2 ( z ) = a/2z/ ny n (z) y 0 (z) = 


-cos z 


Bessel functions 


pfrtir op)- r - f c -u 2 % 2 ( ug cos ^ r ~ + 0 + k a <y) sin M r - gjXgg cos[M(r 0 - g-)] + (i + ^q-)] sin[ M (r 0 - a)] du Substituting 

Jr^l m(l + 2k a <j + (k 2 +u 2 )a 2 ) 


p(r,t|r„,0,0) = 


r 0 V 4;zDt 


exp 


( r ~ r o ) 2 

4Dt 


+ exp 


( r + r 0 ~ 2(7 ) 2 
4Dt 


■ + «W 


V 


r + r 0 - 2cr 

V4Dt 


-6^ 


Vot 


= 4 ^ 2 Prad( r T|r 0 ) 


= Radial GF! 
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Green’s functions (3D) for radiation chemistry 


The sampling algorithms for chemical reactions yields values 
of r, 0 and <|> 

Therefore, we need to work with the distributions p(r,t\r 0 ,6 0 ,</> 0 ), 
p(6,t\r o ,6 o ,0 o ) and p(<p,t\r 0 ,6 0 ,<p 0 ). 

These distributions are difficult to obtain for the general case 
(0 o *O, (j^O) even for simple free 3D diffusion 


Inn 


Pfr(r.t|r o ,0 o .A)= J|r 2 Pfr(r. 0 , 4 , 1 1 r o ,0 o ,$,)sinftifti0 = -j==exp 


r 2 +r 2 ^ 


0 0 

oo In 




4Dt 


sinh 


^ rr ft A 


v 2Dty 


p fr ( 6 »,t | r 0 O 0 ,<j> 0 ) = J jVp fr (r, 6 ^,t| r o ,0 o ,0 o )sin(0)d fair 


0 0 


( 2 w H~ 2 )!^ _^( 1 _ C 0 S ^ C 0 S ^ 0 o )/ 4 Dt J - 2 (n+l) 


n = 0 

CO K 


w 


i Mn+,) Erfc 


- r n cos 0 cos 0, 


V 


r n sin 0 sin O 0 
y[4Dt 


\2 n 


sin# 


P fr 1 1 r 0 O 0 , <f > 0 ) = f f r 2 p fr (r, 0, 1 1 r 0 , # 0 , ^ )sm(0)d0dr = —= 
oo 2 V/r 


V4Dr 

1 - 0 2 sfa 2 (<Mo)sm 2 ^/ 4 fl f f r 0 cos(^ - ^ 0 ) sin # 0 


iErfc 


jADt 


y 


GO 

Erfc(z) = | i n ~ l Erfc(t)dt 


z 
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Green’s functions (3D) for radiation chemistry 


• The sampling algorithms for chemical reactions yields 
distributions of r, 0 and § 

• The analytical distributions p(r,t\r 0 ,6 0 ,<f> 0 ), p(G,t\r 0 ,Q 0 ,<f> 0 ) and 
p(<f>,t\r 0 ,Q 0 ,</> 0 ) are difficult to obtain for the general case (0 o *O, 
<j) 0 ^0) even for simple free diffusion 

• The case 0 o =O, (|) 0 =0 is much simpler since <(> is uniformly 
distributed between 0 and 2n. 

• It is possible to sample the distributions in a referential aligned 
with the z axis and then perform two rotations to get the 
sampled values for the general case 


X 


cos^ 0 

- sin (/) Q 

0 " 

cos# 0 

0 

sin# 0 

x" 

y 

— 

sin^ 0 

cos^ 0 

0 

0 

1 

0 

/' 

z 


0 

0 

1 

- sin 0 0 

0 

COS0 o 

z" 
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Green’s functions (3D) for radiation chemistry 


Free diffusion for the case 0 o =O, 4> 0 =0 

The distributions in Cartesian coordinates are Gaussians. 
Their equivalents in spherical coordinates are given by 


p fr (r,t|r 0 ,0,0) 


Pfr(^t|r 0 ,0,0) 


V 


;zr 0 2 Dt 


exp 


r +r 0 
4Dt 


2 A 


sinh 


rro 

2Dt 


1 

2k 


Pfr (#> 1 1 r 0 ,0,0) = 2(sin 6)e 


-c 0 2 sin 2 6! 4Dt 





r 0 cos 6 ^ 

2 -\[l)t y 


Cartesian coordinates 


Spherical coordinates 



Coordinate 
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Green’s functions (3D) for radiation chemistry 


We know how to sample p(r,t\r 0 ,6 0 ,<f> 0 ). 

The exact Green’s function is difficult to evaluate numerically. 

To calculate it we write p(r,6),<f>,t\r 0 ,6 0 ,</> 0 ) as the sum of free 
diffusion and correction terms (van Zon and ten Wolde, 2005) 

This approach reduces the number of terms necessary for 
convergence 


Pex ( r >0><M I r 0 , 0 O , </> 0 ) = p fr (r, 0, <f>, 1 1 r 0 , 0 O , ip 0 ) + p C0II (r, 6, </>, 1 1 r 0 , 0 O , </> 0 ) 


I yi 00 

'nnax - 

Pcorr ( r > ^ ^ 1 1 r 0 , 0 O , <p 0 ) = - - — -7= £ (2n + l)P n (cos 0) J e‘ u 


Dt 


R, 


4/rJn u „, 0 


+ R\ 


(R,G, + R 2 G 2 )udu 


R, =(2 ak a + \)J, Mjl (ua)-2uaJ', Mjl {ua) 


^1 ^ \ i +\!2 Y n + 1/2 (“'') 7 » + 1 / 2 ( W '' o ) 


R 2 =( 2ak a +l)Y„ +in (u(j)-2uo-Y’ ll+l/2 (u(j) 


^ 2 *^n+l/2 (ur)Y Mll2 (ur 0 ) + 

Y n + 1/2 ( Ur V,Mll( Ur o) 
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Green’s functions (3D) for radiation chemistry 


To sample 0 and <|>, we generate values for 0 o =O, (|> 0 =0 and perform the 
necessary rotations. In this system, <|> is uniformly distributed between 0 
and 2ji 

Use the algorithm from Clifford et al. to generate the angle 0 and (|> 



1.2 


1 . 0 - 



where U 1 and U 2 are random 
numbers uniformly distributed [0,1] 
a=rr n /2Dt 


m — I I q/ Ly l 

r 0 and r are the inter-particle 
distances before and after the 
time step 


0 . 8 - 


0 . 6 - 


0 . 4 - 


0 . 2 - 


0.0 


0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

0 (rad) 
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Radiation chemistry 


T«*r(4ffi5e» 
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Radiation chemistry 
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G x (Molec./lOO eV) 


Radiation chemistry 


Primary yields of e - aq> OH, H , H 2 and H 2 0 2 as a function of the LET 



Irradiation by 300-0.1 MeV protons 
LET: -0.3-85 keV/jum 


Note: the / primary, yields (noted G x ) are 
the yields at the end of spur expansion 
(~10' 6 s) 
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Radiation chemistry (DNA) 


□ The radiation chemistry of DNA is very complex 

□ Many reaction rate constants are known 



Reaction 

k 

Radius 


(dm 3 .mol -1 .s _1 ) 

(A) 

e "aq + Thymine— >Thy(+e) 

1.79x10 10 

5.287 

■OH + Thymine —*■ 
TC50H+TC60H+TUCH2- 

6.4x1 0 9 

3.02 

H.+Thymine— >Thymine* 

5.7x1 0 8 

0.11 
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Release history 



□ RITRACKS was used by the 
students at the NASA Space 
Radiation Summer School at the 
Brookhaven National Laboratory, 
Upton, New York (June 6-24, 2011, 
May 28 - June 15, 2012, MIT ICED 
June 2012)(over 40 users) 

□ The release to international partners 
was approved in 2011 

□ RITRACKS was released to NASA 
space radiation community with 
over 20 users 

□ The software is now available for 
download on the web site 
http://spaceradiation.usra.edu/irMod 
els/ (ITAR, authentication and 
password required) 



http://spaceradiation.usra.edu 
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Release history 


□ An online version of RITRACKS will be available soon! 


Ion; 12C6+, 60 NeV/anu Track No;G Ion; 12C6+, 60*0000 NeV/anu Track No;0 




Left: Simulation of a 12 C 6+ , 60 MeV/amu, on the projected RRAW site. 

Right: Calculation of the radial dose for the track depicted on the left 

The track structure data and the radial dose are available for download after calculation 
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Future plans for development and use 


□ Implementation of the non-homogeneous chemistry 

□ Predictions of clustered and complex DNA damage 
yields in human cells for improving the 
understanding of DNA repair and signal 
transduction 

□ Use with chromosome models to study double- 
strand breaks (DSB) in relation to cancer risks from 
space radiation 

□ Web-based version 

□ New GPU-CPU version to improve computational 
speeds by several orders of magnitude. 
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Future plans for development and use 



□ Android/iPhone version 

□ For ions: 

• LET 

• Relativistic p and y 

• Z eff and z eff 2 /p 2 

• Maximum energy transfer to an electron 

• Dose and fluence 

• Radial dose 

• Number of hits per cell in a cell culture 

□ For electrons: 

• Relativistic (3 and y 

• Range 
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Future plans for development and use 



□ Android/iPhone version 

□ For ions: 

• LET 

• Relativistic (3 and y 

• Z eff and Z eff 2 /p 2 

• Maximum energy transfer to 
an electron 

• Dose and fluence 

• Radial dose 

• Number of hits per cell 

• in a cell culture 

□ For electrons: 

• Relativistic (3 and y 

• Range 
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Using RITRACKS 



□ The installer 

• The necessary files are included in an installer created by the freeware 
InstallJammer 

• The program is installed in the folder C:\Program Files (x86)\RITRACKS 

• Simulations are stored the subfolder RITRACKS Simulations in the My 
Documents folder 



RITRACKS Setup 


Choose Destination Location 

Where should RITRACKS be installed? 



Setup will install RITRACKS in the following folder. 

To install to this f older , click Next. To install to a different f older , click Browse and select another 
folder. 


Destination Folder 
C:\Prograrm Files (x86)\RITRACK5 


Browse... 


InstallJammer - 


Next > 

| | Cancel 


< Back 

Next > ;| 


Cancel 
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RITRACKS main window 
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Radiation info window 



The following information is given in this window 


Rest mass energy: 

Total energy: 

Relativistic y: 

Relativistic p: 

Momentum: 

Maximum energy 
Transfer to e - : 


Me 2 

yMc 2 

y = T/Mc 2 +1 
P 2 = 1-1/y 2 


p = yMv 


2mc 2 (y 2 -1) 

1 + 2y(m / M) + (m / M) 2 



LET (MeV/cm): 
(Bethe) 


dE 

dx 


0.170 

~P^ 


[p(P) - 4.3 1] 


F(P) = In 


1 .02x1 0 6 p 2 

1-P 2 


P 2 
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Cross sections 


Electron cross sections window 




Number of samples: 


10000 
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Visualization window 



Tools: - Rotation 

- Translation 

- Zoom 

- Save to file 

- Copy to clipboard 

- Create a .avi file 

- Open data folder 

Visualization: 

- Radiolytic species 

- Events 

- Dose (voxels) 

Time evolution 
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RITRACKS tools 



□ Calculation of tracks per cell in a cell culture for a 
given ion, energy and dose 
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Homework 



□ Install RITRACKS 

□ Max the number of CPUs 

□ Click on the More button for the 
Radiation Info window 



Use this form to select the CPU ressources used for the calculations 

Number of NUMA nodes: 2 

Number of physical processor packages: 2 

Number of processor cores: 12 

Number of logical processors: 24 

Number of CPU used for the calculation: 24 

Processor details 
Use No 



Close 
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Homework 



□ Simulations: 

• 1 ion 12 C 6+ , 25 MeV/amu, of length 5 jum. Use the following 
parameters: 

S #Copy = 1 and #Histories/Copy = 20 

S #Copy =10 and #Histories/Copy = 2 

S Notice the difference in the progress window and calculation time 

□ Visualize the tracks with the 3D interface 

• Use visualization tools (displace, rotate, zoom, copy to clipboard, save 
to file, reset parameters, open folder). 

• Zoom and un-zoom (to un-zoom, draw a zoom box from the bottom 
right to upper left) 

• Use the time evolution bar. 

• Change the simulation number and click Draw. 

• Create an avi file of the track. Click on the movie icon, click on Create 
AVI, enter a name (test.avi). 
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Homework 



□ Visualize the tracks with the 3D interface 

• Switch between voxel dose and track. In the voxel dose 
representation, vary the threshold scale (the slider located at the left of 
the dose scale) 

• Display the axis (go to Options->Visualization, select the tab Axis and 
check the boxes). 

□ Display the radial dose and voxel dose (from the main 
interface window). 

□ Repeat the simulations with the following ions 

• 1 ion 56 Fe 26+ , 1000 MeV/amu, of length 5 pm (10 histories) 

• Irradiation of a cube of 5 pm x 5 pm x 5 pm with 450 x 300 MeV 
protons (1 history) 

• Irradiation of a cube of 5 pm x 5 pm x 5 pm with 1 1000 MeV/amu 
56pe 26+ ions (1 history) 
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